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Foreword 


The  Naval  Oceanographic  and  Atmospheric  Research  Laboratory’s 
Mapping,  Charting,  and  Geodesy  Division  performs  research  and  develop¬ 
ment  in  new  and  improved  methods  for  gathering  navigational  bathymetric 
data  in  coastal  areas.  One  area  of  special  interest  is  the  development  of 
algorithms  to  derive  bathymetry  from  multispectral  data.  Currently,  multi- 
spectral  imagery  can  be  derived  from  airborne  scanners  or  from  satellite  sensors, 
such  as  the  thematic  mapper  on  board  Landsat.  Multispectral  bathymetry  has 
the  potential  for  gathering  data  much  more  rapidly  than  present  shipboard 
systems;  in  particular,  satellite  systems  can  gather  data  from  many  parts  of 
the  world  that  are  denied  to  other  collection  systems.  This  report  presents 
results  in  the  use  of  different  sources  of  multispectral  data  for  bathymetry 
derivations. 
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Executive  Summary 


The  Mapping,  Charting,  and  Geodesy  Division  of  the  Naval  Oceanographic 
and  Atmospheric  Research  Laboratory  (NOARL)  investigates  methods  to 
exploit  the  application  of  multispectral  data  to  coastal  water  mapping  and 
charting  requirements.  The  Coastal  Image  Understanding  project  high- 
resolution,  multispectral,  aircraft  data  has  been  examined  and  compared  with 
lower-resolution  satellite  data  from  the  thematic  mapper  on  board  NASA’s 
Landsat.  The  high-resolution,  nine-band  scanner  used  was  developed  for  the 
Airborne  Bathymetric  Survey  System.  The  NOARL  scanner  has  increased 
spectral  sensitivity  in  the  blue  and  green  portions  of  the  electromagnetic 
spectrum  to  take  advantage  of  the  water-penetrating  bands.  The  multiband 
generalized  ratio  model  is  used  to  compare  bathymetry  estimates  from  the 
scanner  and  the  thematic  mapper,  which  has  five  bands.  Results  show 
improvement  from  increased  resolution  in  the  3-m  to  6-m  area.  Root-mean- 
square  (RMS)  errors  of  0.323  m  at  the  3-m  depth  range  for  the  NOARL  Scanner 
are  close  to  the  International  Hydrographic  Office  standard  of  0.3  m;  however, 
errors  at  other  depth  ranges  exceed  this  standard.  RMS  errors  of  0.544  m  for 
the  NOARL  scanner  contrast  with  0.907  m  for  the  thematic  mapper  at  the 
4-m  range.  Beyond  6-m  depth,  results  do  not  differ  significantly  between 
the  two  scanners. 
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Comparison  of  Thematic  Mapper  and 
NOARL  Scanner  Multispectral  Bathymetry  Estimates 


Introduction 

Significant  advances  have  been  made  in  the  use 
of  satellite  multispectral  imagery  for  the  estimation  of 
bathymetry  (Clark  et  al.,  1988;  Clark  et  al.,  1987; 
Lyzenga,  1985).  Extensive  analysis  of  thematic 
mapper  (TM)  imagery,  gathered  on  board  NASA’s 
Landsat,  has  shown  the  multiband  regression  model 
to  be  robust  and  accurate  to  within  1  to  2  m.  Develop¬ 
ment  of  the  airborne  NOARL  scanner  provided 
multispectral  data  with  both  higher  ground  and  spectral 
resolution.  The  TM  has  a  30-m  ground  resolution, 
whereas  the  scanner  has  a  1-m  ground  resolution  (based 
on  altitude  of  500  m).  Table  1  compares  the  scanner 
and  TM  spectral  bands.  Increased  resolution  in  TM 
bands  1  and  2  provides  more  information  for  the 
scanner  in  the  water-penetrating,  blue-green  bands. 

A  major  goal  of  the  Coastal  Image  Understanding 
project  is  to  compare  the  effects  of  the  scanner’s 
increased  resolution  on  bathymetry  estimates  with  those 
obtained  from  TM  imagery.  This  report  documents  the 
bathymetry  results  using  multispectral  data  from  each 
of  the  sensors  compared  to  ground  truth  data  obtained 
from  a  Naval  Oceanographic  Office  (NAVOCEANO) 
ship  survey  off  Key  West,  FL  (Fig.  1). 


Background 

•  The  Coastal  Image  Understanding  project  is  an 

effort  to  improve  the  accuracy  of  passive  shallow-water 
bathymetry  and  to  classify  and  identify  features  of 


Table  1.  Spectral  bandwidths  of  the 
scanner  and  the  TM. 


SCANNER 

TM 

BANDWIDTH 

1 

1 

450 

-  480  nm 

2 

1 

480 

-  520  nm 

3 

520 

-  550  nm 

4 

550 

-  600  nm 

5 

630 

-  690  nm 

6 

760 

-  900  nm 

7 

5 

1.55 

-  1.75  fim 

6 

6 

2.08 

-  2.35  \k m 

9 

7 

10.40 

-  12.50  fim 

naval  interest  in  the  coastal  area  that  can  directly 
address  U.S.  Navy  and  Defense  Mapping  Agency 
(DMA)  requirements. 

Past  work  on  TM  data  (Clark  et  al.,  1988)  has  shown 
that  the  multiband  algorithm  performs  better  than  the 
previously  tested  single-  and  dual-band  algorithms. 
The  multiband  regression  model  using  TM  gives  1-  to 
2-m  root-mean-square  (RMS)  depth  errors  in  the 
relatively  clear  waters  off  Key  West,  FL,  and  the  island 
of  Puerto  Rico. 

The  development  of  the  high-resolution  multispectral 
scanner  for  the  Airborne  Bathymetric  Survey  (ABS) 
system  produced  data  with  1-m  resolution  and  finer 
bandwidths  in  the  water-penetrating  portion  of  the 
visible  spectrum.  During  the  field  testing  of  the  ABS 
system,  much  data  was  collected  over  the  Key  West 
area  around  the  East  Sambo  Shoals,  including  dense 
hydrographic  boat  soundings,  as  well  as  multispectral 
data  over  the  same  area. 

The  application  of  the  multiband  algorithm  to  the 
high-resolution  scanner  data  is  of  interest  in  deter¬ 
mining  whether  that  data  can  give  better  depth 
estimates.  The  availability  of  scanner  data,  along  with 
TM  data  and  bathymetry  data  in  the  Key  West  area, 
provided  the  data  set  for  comparing  the  various  resolu¬ 
tion  effects  on  depth  estimates. 


Procedures 

Current  multispectral  depth  algorithms  are  based  on 
the  functional  form  of  the  single-band  radiance  model 
(Clark  et  al.,  1987a).  Statistical  regression  techniques 
are  used  to  determine  the  unknown  parameters  using 
water  depths  known  at  several  pixel  locations  in  the 
multispectral  image.  Single-  and  dual-band  ratio 
models  are  extremely  sensitive  to  changes  in  bottom 
reflectance  and  water  clarity.  Multiband  algorithms  in 
which  the  radiance  and  bottom  reflectance  are 
functions  of  the  wavelength  can  be  used  to  minimize 
difficulties  encountered  due  to  changes  in  the  bottom 
reflectance  and  water  attenuation.  A  development  of 
the  multiband  model  from  the  single-band  radiance 
equation  is  presented  here. 


1 


2 


The  single-band  radiance  equation,  where  total 
radiance  from  the  water,  Lt,  for  a  single  bandwidth 
(/)  is  written  as 

Lt  =  Lmj  +  S'  rt  expl-fkjZ],  (1) 

where 

Leo i  =  radiance  from  deep  water; 

S'  =  a  function  of  transmittance  of  air  and  water 
surface,  surface  refraction,  and  solar  irradiance; 

r(  =  bottom  reflectance; 

k )  =  effective  attenuation  coefficient  of  the  water; 

/  =  geometric  factor  to  account  for  the  pathlength 
through  the  water; 

z  =  water  depth; 

i  =  1,2  ...  n; 

n  =  number  of  spectral  bandwidths. 

This  model  has  had  several  implementations  (Clark 
et  al.,  1987b;  Paredes  and  Spero,  1983).  Some 
implementations  require  either  classifying  the  data  into 
bottom  cover  types  or  knowledge  of  optical 
parameters.  Other  implementations  use  band  ratios  and 
combinations  of  bands.  Ratios  work  rather  well  if  they 
are  constant  for  all  bottom  types  in  the  area  of  interest. 
The  ratio  method  was  outlined  by  Polcyn  et  al.  (1976) 
and  elaborated  by  Lyzenga  (1985).  The  generalized 
band  ratio  algorithm  has  been  found  to  give  the  best 
results  (Clark  et  al.,  1988).  A  brief  summary  of  the 
method  is  given  here.  Equation  (1)  is  usually  linearized 
and  solved  for  z  to  yield 

Z  =  (l//k,)(/«[r,]  -X,),  (2) 

where  X,  =  (/j[(L,  -  Loo, )/$,].  Equation  (2)  can 
be  reparameterized  with  bai  =  (1  /fkf)  (/«[/•,]), 
bu  =  -(1 //*,),  and  X0i  =  1 

The  depth  z,  for  1  channel,  can  the  be  expressed  in 
matrix  notation  as  z  =  b  'x.  The  vector  b  contains  the 
weights  bQ  and  br  The  vector  x  contains  [1,  x]. 

For  n  bandwidths,  the  depth  z  can  be  expressed  as 
a  weighted  sum  of  single  band  terms  as  follows: 

z  =  ^,6,*,  +  w2b2x2  +  ...  +  wn  b„xn,  (3) 

where  the  b  and  x  vectors  are  defined  as  above  for  each 
band,  and  the  ws  are  weights  whose  sum  equals  1.  In 
matrix  form,  this  equation  can  be  written 

Z  =  w'B'x,  (4) 

where  w  is  a  vector  of  n  +  1  weights,  B  is  a  diagonal 
matrix  of  weights  where  b0  is  a  constant  term  derived 
from  the  sum  of  constant  terms  for  each  band,  and 
x  is  an  n  +  1  vector  with  x0  =  1.  The  motivation 
behind  the  summation  of  constant  terms  is  given  by 
Paredes  and  Spero  (1983),  and  discussed  here. 


For  m  bottom  cover  types,  and  n  bandwidths  such 
that  m  =  n,  Paredes  and  Spero  have  shown  that  a 
vector  c  of  length  n  exists,  such  that 

c'A  =  1,  (5) 

where  1  is  a  1 '  s  vector  of  length  n  and  A  is  defined 
as  the  n  by  m  matrix. 

In  rn  ....  In  rlm 

A  =  In  r2l  ....  In  r2m  (6) 


Each  element  of  A  is  the  logarithm  of  the  bottom 
reflectance  term  for  the  nth  channel  and  the  mth 
bottom  type.  As  long  as  the  number  of  bottom  types, 
m,  is  not  greater  than  the  number  of  channels,  n,  a 
solution  for  c  can  be  found.  Paredes  and  Spero  show 
that  for  m  =  n  bottom  types,  the  weights  w,  can  be 
written 

vv(  =  Cjkj  /  c'k.  (7) 

Equation  4,  can  be  written  as 

z  -  w’b0  +  w'bix,  (8) 

where  w'b0  =  (l//)[c'lc  c’A.  Since  c'A  is  1,  the 
bottom  reflectance  term  vanishes,  and  only  the  sum 
of  constant  terms  remains.  The  multiband  depth  model 
can  thus  be  expressed  independently  of  bottom  type 
reflectance,  as  the  linear  regression  model 

z  =  fiX  +£,  (9) 

where  fi  is  the  vector  w'B  of  coefficients  derived  from 
the  least-squares  solution  of  equation  3,  X  is  the  vector 
of  reflectances,  and  £  is  the  error. 

Results  from  comparison  of  the  multiband  model 
to  the  single-band  and  dual-band  ratio  models  indicate 
that  the  multiband  generalized  ratio  technique  gives 
superior  results  and  is  computationally  efficient  as  well. 

Results 

Comparisons  were  made  using  various  channel 
combinations.  In  comparing  these  different  spectral 
bandwidths  for  bathymetry  estimation,  it  is  useful  to 
observe  the  transmittance  of  light  through  pure 
seawater  (Fig.  2).  The  transmittance  falls  off  very 
quickly  after  575  nanometers  (nm).  The  transmittance 
at  wavelengths  below  575  nm  is  greater  than  90%  for 
pure  seawater.  The  wavelength  range  with  the  higher 
transmissivity  will  therefore  be  the  best  for  depth 
determination.  Transmittance  for  coastal  waters  shows 
a  similar  dependence  but  is  of  lower  magnitude.  The 


range  of  wavelength  bands  for  the  two  scanners  is 
depicted  by  vertical  bars  in  Figure  2. 

Tests  were  also  run  using  a  maximum  depth  range 
to  examine  the  depth  effects  on  the  bathymetric  error 
estimates  from  the  different  scanners.  Different  band 
comparison  results  are  presented  as  a  function  of 
maximum  depth.  The  depth  range  used  for  each 
comparison  always  ranges  from  0  m  to  maximum 
depth. 

The  test  data  set  includes  ship-derived  depth 
soundings  from  a  hydrographic  survey  by 
NAVOCEANO  conducted  as  part  of  the  ABS  system 
field  tests  in  August  of  1988.  The  survey  area  extended 
from  near  24°29.5'  to  24°31'N  and  from  81°39.75' 
to  81 0 39.25 '  W.  The  area  (Fig.  1)  covers  the  eastern  half 
of  the  Eastern  Sambo  shoal  area  and  the  region  to  the 
north  of  it.  Approximately  540  depth  soundings  exist 
for  the  scanner  and  TM  data  sets.  The  depths  varied 
from  0  to  10  m.  The  spacing  between  boat  soundings 
(25  m)  was  dense  due  to  the  high  ABS  sampling  rate. 
The  TM  data  in  most  cases  had  more  than  one 
sounding  per  30-m  pixel.  In  these  cases,  the  soundings 
were  averaged  to  produce  a  mean  depth  per  pixel. 

Figure  3  shows  the  results  of  the  multiband  model 
using  all  visible  bands  of  the  TM  and  all  visible  bands 
of  the  scanner.  Table  2  summarizes  the  results.  The 
errors  were  slightly  less  for  the  scanner  in  the  3-  to  6-m 
range.  The  differences  extend  from  0.36  m  for  waters 
up  to  4  m  deep  to  0.03  m  for  waters  up  to 
6  m  deep.  Differences  for  depths  beyond  6  m  were  not 
significant.  Both  data  sets  showed  an  upward  trend 


350  400  450  500  550  600  650  700  750  800 
WAVELENGTH  (nm) 

Figure  2.  Optical  transmittance  of  pure  seawater  versus 
wavelength  (after  Jerlov,  1968).  TM  and  scanner  bandwidths 
shown  by  vertical  lines.  (Transmittance  for  coastal  waters 
shows  a  similar  dependence  but  is  between  5  and  10%  lower 
in  magnitude .) 


MAXIMUM  CALIBRATION  DEPTH  (m) 

Figure  3.  Comparison  of  errors  in  bathymetry  estimation 
between  scanner  bands  1,  2,  3,  4,  5  and  TM  bands  1,  2,  3. 


Table  2.  Summary  of  scanner  (S)  and  TM  bathymetry  estimates  using  TM  bands 
1,  2,  3,  and  S  bands  1,  2,  3,  4,  5.  Errors  in  meters. 


in  the  RMS  error  with  increasing  depth.  Scanner  error 
increased  steadily  and  leveled  off  around  9  m.  TM 
error  seemed  to  increase  more  rapidly  in  the  shallower 
waters  between  3  and  5  m  and  followed  the  same  trend 
as  the  scanner  after  6  m. 

The  mean  residual  error  is  plotted  in  Figure  4.  The 
TM  showed  a  consistent  negative  bias.  Since  the  errors 
are  cumulative,  it  was  believed  the  errors  at  the  3-m 
level  were  weighting  the  errors  at  the  higher  depths 
downward;  however,  the  error  at  3  m  was  based  on 
only  6  points.  The  error  at  10  m  is  based  on  147  test 
points.  It  is  possible  that  the  La,  values  used  were  too 
low.  The  errors  for  the  scanner  were  less  erratic  and 
less  biased  and  were  constrained  between  0  and  0.08  m. 
The  errors  for  the  TM,  however,  were  between  -0.05 
and  -0.20  m  (excludes  3-m  error).  The  errors  did  not 
show  trends. 

In  addition  to  the  all-band  comparisons,  different 
band  combinations  were  compared.  As  shown  in 
Table  1 ,  scanner  bands  1  and  2  cover  the  same  spectral 
bandwidth  as  TM  band  1 ,  and  scanner  bands  3  and  4 
cover  the  same  spectral  bandwidth  as  TM  band  2. 
Scanner  band  5  is  equivalent  to  TM  band  3;  therefore, 
combinations  of  the  scanner  bands  were  tested  against 
the  TM  bands.  Test  results  are  shown  in  Tables  3  and  4, 
and  error  plots  are  shown  in  Figures  5  and  6.  Results 
of  the  two-band  and  four-band  runs  for  the  scanner 
did  not  differ  greatly,  mainly  due  to  the  high  correla¬ 
tion  between  the  blue  and  green  bands.  Figure  7  shows 
scatter  plots  of  various  band  combinations  in  which 
the  high  degree  of  correlation  can  be  seen.  The  best 
results  from  the  two-band  combinations  came  from 
bands  2  and  4.  The  five-band  results  were  slightly  but 
not  significantly  better  than  the  two-band  results.  Com¬ 
parison  with  the  TM  bands  showed  similar  results  to 
the  all-band  comparisons,  with  better  results  from  the 
scanner  in  the  3-  to  6-m  range,  and  no  significant 
differences  beyond  7  m. 

Figure  8a  shows  the  TM  true-color  image  with  the 
scanner  coverage  area  outlined  in  red  (Fig.  8b). 
Figure  9  shows  the  true-color  image  of  the  raw  scanner 


Figure  4.  Plot  of  mean  residual  error  for  scanner  and  TM. 


data  for  the  area  outlined  in  Figure  8b.  Figures  10 
and  1 1  show  bathymetry  contours  from  the  TM  and 
scanner  data  sets,  respectively. 

The  wider  scanner  bands  that  correspond  to  the  TM 
are  used  to  compare  ground  resolution  effects.  Scanner 
bands  2  and  4  are  compared  with  TM  bands  1  and  2. 
The  result  is  the  same  as  that  reported  for  the  two- 
band  comparisons  shown  in  Figure  5. 

Summary  and  Conclusions 

TM  and  scanner  data  were  compared.  The  scanner 
data  were  better  than  the  TM  data  in  the  shallow  area 
between  3  and  6  m.  Comparisons  of  the  scanner  data 
at  different  resolutions  is  planned  to  differentiate  the 
combined  effect  of  increased  spectral  and  spatial 
resolution. 

The  improvement  in  data  by  using  the  scanner  was 
less  than  expected  and  is  attributed  to  several  factors. 
First,  the  scanner’s  higher  spatial  resolution  is  much 
more  sensitive  to  horizontal  registration  errors 
than  the  lower  resolution  TM.  The  integrated  response 
over  the  lower  resolution  TM  and  the  averaging  of 
soundings  within  the  30-m  cell  area  could  give 


Figure  5.  Comparison  of  errors  in  bathymetry  estimation 
between  scanner  bands  2,  4,  and  TM  bands  I,  2. 


Figure  6.  Comparison  of  errors  in  bathymetry  estimation 
between  scanner  bands  4,  5  versus  TM  bands  2.  ?. 


Table  3.  Summary  of  scanner  (S)  and  TM  bathymetry  estimates  using  TM  bands  1  and  2  and  S  bands 
1,  2,  3,  and  4.  Errors  in  meters. 
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-0.4167E-01 

0.3200 

12 

12 

-0.9537E-06 

0.3227 

12 

12 

-0.8333E-01 

0.3118 

12 

12 

-0.9537E-06 

0.3227 

12 

12 

-0.1167 

0.3436 

6 

6 

-0.9537E-06 

0.5701 

38 

40 

0.8750E-01 

0.5743 

38 

40 

-0.1250E-01 

0.4998 

38 

40 

0.5000E-01 

0.5679 

38 

40 

-0.1167 

0.6700 

14 
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0.1027E-01 

0.7590 

73 

73 

0.7877E-01 

0.6785 

73 

73 

-0.1027E-01 

0.8516 

73 

73 

0.7877E-01 

0.7225 
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73 

-0.1923 

0.8358 

24 

26 

-0.6429E-01 

0.9790 
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-0.3929E-01 

0.8368 

142 

140 

-0.2143E-01 

1.204 
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-0.1429E-01 

0.8833 
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-0.6035 E-01 
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93 
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234 

0.6624E-01 

1.186 

236 

234 

0.1068E-01 

1.345 

236 

234 

0.5555E-01 

1.159 

236 

234 
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1.120 

118 
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0.6509E-01 

1.459 

265 

265 

0.6887E-01 

1.293 

265 

265 

0.4245E-01 

1.485 

265 

265 

0.3113E-01 

1.243 

265 

265 
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1.245 

143 
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0.6227E-01 

1.484 

270 
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1.301 

270 

269 

0.4740E-01 

1.514 

270 

269 

0.6041  E-01 

1.269 

270 

269 

-0.1378 

1.285 

147 

147 

improved  estimates  on  the  order  of  the  square  root  of 
the  number  of  soundings.  The  higher  number  of  bands 
does  not  seem  to  offer  much  improvement  due  to  the 
high  amount  of  collinearity  between  them.  However, 


the  finer  spectral  bandwidths  can  be  used  to  advantage 
in  different  types  of  ocean  water,  since  the  attenua¬ 
tion  and  absorption  properties  are  optimal  at  different 
wavelengths. 


6 


Table  4.  Summary  of  scanner  (S)  and  TM  bathymetry  estimates  using  TM  bands  2  and  3  and  S  bands 
3,  4  and  5.  Errors  in  meters. 


DEPTH 

SENSOR/BANDS 

RMUL 

{ACTUAL-CALCULATED) 
MEAN  RMS 

#  CALIB 
POINTS 

#  TEST 
POINTS 

0-3 

S/35 

0.574 

-0.9537E-06 

0.3227 

12 

12 

S/45 

0.756 

-0.4167E-01 

0.3200 

12 

12 

TM/23 

0.931 

-0.3333 

0.4488 

6 

6 

0-4 

S/35 

0.165 

-0.9537E-06 

0.5244 

38 

40 

S/45 

0.305 

-0.1250E-01 

0.5122 

38 

40 

TM/23 

0.614 

0.1167 

0.8260 

14 

15 

0-5 

S/35 

0.207 

-0.3767E-01 

0.7935 

73 

73 

S/45 

0.179 

-0.2397E-01 

0.8111 

73 

73 

TM/23 

0.359 

-0.1154 

0.9152 

24 

26 

0-6 

S/35 

0.211 

-0.6786E-01 

1.504 

142 

140 

S/45 

0.311 

-0.1179E-00 

1.006 

142 

140 

TM/23 

0.617 

-0.1638 

1.153 

58 

58 

0-7 

S/35 

0.273 

-0.7896E-02 

1.217 

190 

190 

S/45 

0.357 

-0.2632E-01 

1.153 

190 

190 

TM/23 

0.550 

-0.2957E-01 

1.177 

82 

93 

0-8 

S/35 

0.322 

-0.1710E-01 

1.424 

236 

234 

S/45 

0.430 

-0.2351  E-01 

1.373 

236 

234 

TM/23 

0.410 

-0.1618 

1.358 

118 

119 

0-9 

S/35 

0.285 

0.3868E-01 

1.651 

265 

265 

S/45 

0.362 

-0.9444E-03 

1.584 

265 

265 

TM/23 

0.158 

-0.9722E-01 

1.705 

143 

144 

0-10 

S/35 

0.281 

0.1208E-01 

1.671 

270 

269 

S/45 

0.348 

-0.4648E-02 

1.621 

270 

269 

TM/23 

0.142 

-0.1173 

1.698 

147 

147 

CHANNEL  3  VS.  CHANNEL  1 


Figure  7.  Scatter  plots  of  scanner  channels  showing  the  correlation  between  the  blue  and  green  bands  for  different 
combinations. 


Figure  8.  (a)  TM  image  of  Key  West,  FI.  (b)  Study  area  outlined  in  red. 


-e  9.  True-color  image  of  Key  West,  FL,  from  scanner. 
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Recommendations 

Recommendations  for  future  work  are  to  collect  data 
in  an  area  with  a  homogeneous  bottom,  such  as  sand, 
to  minimize  errors  due  to  rapidly  changing  bottom 
types.  The  registration  of  1-m  resolution  scanner  data 
to  a  rapidly  changing  bottom  requires  very  high  preci¬ 
sion  positioning  data.  The  limited  availability  of  GPS 
satellites  during  the  collection  of  the  scanner  data  did 
not  yield  the  best  satellite  configuration  for  the  precise 
positioning  needed.  A  slowly  changing  bottom,  where 
most  of  the  variability  is  from  depth  variations  and 
not  from  cover  type  changes,  would  be  more  robust 
to  registration  errors. 

Research  also  needs  to  be  done  to  see  if  surface 
characteristics,  derivable  from  glint  patterns  on  the 
surface  in  the  multispectral  data,  can  be  used  to 
improve  bathymetry  estimation.  Surface  waves  can 
provide  some  information  about  the  nature  of  the 
bottom,  which  might  then  be  used  to  help  distinguish 
depth  variations  from  bottom  type  variations.  This 
distinction  could  provide  better  depth  accuracy  from 
the  multiband  model  in  areas  where  bottom  type 
variability  is  a  problem. 

For  areas  of  the  world  where  surveying  is  denied, 
the  application  of  bathymetry  collection  techniques 
from  multispectral  data  could  yield  products  descriptive 
of  coastal  zone  variability,  surf-induced  beach  erosion, 
shoal  emergence  areas,  and  possible  shipping  routes. 
Although  bathymetry  can  be  accurately  described  only 
to  3-m  depths  (by  International  Hydrographic  Office 


standards)  with  the  scanner,  satellite-derived 
bathymetry  can  be  valuable  for  planning  purposes  in 
areas  where  no  data  are  available. 
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